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ABSTRACT 

In this letter, we show that a massive black hole (MBH) which falls into the center of a 
O ! galaxy in dynamical timescale leaves a weak cusp (p oc r~^'^) around it, which is in good 

agreement with the recent observations of luminous ellipticals by Hubble Space Telescope. 
Such event is a natural outcome of merging of two galaxies which have central MBHs. 
This is the only known mechanism to form weak cusps in luminous ellipticals. Therefore, 
the existence of the weak cusps indicates the central BHs of luminous ellipticals have 






i-{ • fallen to the center from outside, most likely during a major merger event. 

a 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: kinematics and 
v~j ■ dynamics — galaxies: nuclei — galaxies: structure 



Introduction 

Recent observations of el liptical galaxies by Hub- 
ble Space Telescope {HST) { Lauer et al. 1995 ;_ Byun 
3t al. 1996|; [Gebhardt et al. 1996|; |Faber et al. 1997| ; 



Kormendy et al. 1996 ) have brought us some chal- 
lenging problems as well as new knowledges on the 
structure of the central regions of elliptical galaxies. 
First, these observations showed that no elliptical has 
"core", in which the surface brightness or the lumi- 
nosity density profiles would be flat. The central re- 
gions formerly regarded as "cores" proved themselves 
to be power-law cusps in all observed galaxies. Sec- 
ond, the distribution of the slopes of these cusps seem 
to be bimodal, one is the group of "core" galaxies with 
weak cusps (p ex r~",0.4 <, a <^ 0.8) and the other 
is the group of "power-law" galaxies with steep cusps 
{p ex r~", a ~ 2). In addition, the slopes of the cusps 
have correlation with the brightness of the galaxies, 
so that bright galaxies tend to have weak cusps. Car- 
oUo et al. (1997) claimed that the distribution of the 
slope is rather continuous in their sample. However, 
their sample is limited in the range of absolute mag- 
nitude, so it is not clear whether their result is real 
or due to selection. 

The models of cuspy stellar systems previously 
studied are classified into two categories — models 
with and without MBH. However, neither can ex- 
plain the origin of the weak cusps in luminous ellip- 
tical galaxies. It was shown that the cusp around 
BH would have the profile p ex r"^'"* when the evo- 



lution is driven by the thermal relaxati on ( Bahcall & 



Wolf 1976|; Shapiro fc Marchant 197S; Cohn fc Kul 



H-ud 1975 



Marchant fc Shapiro 1979; Marchant fc 



Shapiro 1980 ) and p oc 



^^3jT 



(Young 1980; Merritt fc 



^uinlan 1998) or rather steeper (Quinlan, Hernquist, 



fe Sigurdsson 199J;) when the central BH grows adia- 
batically. Hierarchical clustering in CDM cosmogony 
(Navarro, Frenk, & White 1996; Fukushige & Makino 
[1997 ) or dissipationless collapse (Hozumi, Burkert, & 
[Fujiwara 1999 ) might form relatively shallow cusps, 
but they are still significantly steeper than the ob- 
served weak cusps. 

Makino & Ebisuzaki (1996, hereafter ME96) showed 
that a weak cusp {p oc r^", a ^ 1) is formed through 
the merging of two galaxies which have central BHs. 
They also found that the ratio between the size of 
the weak cusp region and the half-mass radius of the 
merger remnant, Tc/rh, is proportional to the ratio 
between the BH mass and the galaxy mass MBii/Mg- 



This result is in good agreement with observations 
( pebhardt et al. 1996| ). However, they did not dis- 
cuss why the cusp was formed. Quinlan & Hernquist 
(1997) and Makino (1997) showed that the BH bi- 
nary in galactic center ejects many stars as it hardens 
and this process can explain the weak cusps in large 
ellipticals, but they made no quantitative prediction 
about cusp slopes. 



In our previous work (Nakano & Makino 1999, 
hereafter NM99), we investigated the dynamical re- 
action of a galaxy to a BH which falls to the center, 
in order to clarify the formation mechanism of the 
cusp. We found that when the massive BH falls to 
the galaxy center, the stars are heated up by the BH 
and the weak cusp {p oc r~^") is formed (Figure Q[). 
This result is independent of the initial orbital an- 
gular momentum and the mass of BH. Thus, we can 
conclude that when a BH (or BHs) falls from outside 
to the center of a galaxy, a central weak cusp is always 
formed. However, it was not at all clear why the cusp 
is formed. 

In this letter, we present the theoretical explana- 
tion of the formation mechanism of the weak cusp. 

2. Why the Cusp Tends to r ^/^ ? 

In the results of our A^-body simulation, we found 
an important feature of the energy distribution func- 
tion N{£) profiles, shown in Figur e p|. The system of 
units we use is the standard unit ( Heggie fc Mathieu 
1986| ), in which the total mass of a galaxy Mgai = 1, 
the gravitational constant G — 1 and the total energy 
of the galaxy i^tot = —1/4. In this unit, the virial ra- 
dius of the galaxy is scaled to unity and the half- mass 
crossing time is 2-\/2. 

In Fig. g, N{£) is practically unchanged from that 
of the initial King model for the runs in which the 
BH is initially placed off-center. Thus, there is no 
star with £ > 3.3. This is in sharp contrast with 
N{£) for the run in which the BH is initially placed 
on-center (dash-dot-dash line). In this case, N{£) has 
a long tail to £ — > oo. By taking such depletion of 
tightly bound stars into account, we can explain the 
existence of the weak cusp as follows. 

The density p{r) is derived from the distribution 
function f{£) as 

*(-r) 



p{r)=4n f{£W2[^{r)-S]d£, (1) 

Jo 

where ^(r) is the depth of the gravitational poten- 



t = Q (King, W^ = Q]_ 

1, v^^^^O.^v^ - 
^0=1' ^t,o=0-25^K- 




Fig. 1. — Density profiles obtained by TV-body simu- 
lations of BH "fall-in" ( [Nakano fc Makino 1999| ). In 
these runs, the BH mass AfsH = l/24Afgai. BH is 
placed at distance of tq = 1 from the center with tan- 
gential velocity vtfl- Note that wk denotes the initial 
Kepler velocity. In all runs, initial radial velocity of 
BH is zero. 



tial at distance r from the center and £ is the spe- 
cific binding energy. Here we assume the galaxy 
is spherically symmetric and the velocity distribu- 
tion is isotropic. When the galaxy has the cen- 
tral BH, "^{r) near the BH can be approximated as 
^(j') ^ GM-qu/t, and ^(r) diverges when r goes to 
zero. As shown in Fig.H, N{E) of the galaxy after 
BH fell down to the center vanishes at some finite 
value of the binding energy. We denote this limit of 
binding energy as £q. Thus f[£) also vanishes at £o, 
since f{£) — N{£)/A{£) where A{£) is t he area of 
hypersurface in phase space with energy £ ( Binncy & 
Tremaine 19871) . Thus, if *(r) > £o, cq. (yj) can be 



rewritten as 



p{r) = 47r / fi£)^2[-i'{r)-£]d£. 



(2) 



In the region ^(r) > £o, Eq. (0) can be expanded 



as 



£o 



p{r) = An / f{£)^2[^{r)-£]d£ 

JO 



= 4V27rV*W 

X / I{£) 
Jo 



^4^+^ 



*(r) 



d£ 



VW) 



GMbh 



(3) 



Therefore, in the central region where ^(r) ^ £, p{r) 
is proportional to r~^/^. 

Note that our theory is related to, but not the same 
as, the theory for r~^/^ cusp by Zel'dovich & Novikov 
(1971) and Peebles (1972). They showed that if a 
massive object is placed in a stellar system of uni- 
form density (such as very large core), there will be 
small cusp with p oc ^Jl + {GMBH'm/rEoo), where 
Eoo is the energy of a field star. In their theory, this 
cusp is due to gravitational focusing of stars of uni- 
form background. In our theory, the cusp is also due 
to stars which are not bound to BH, but we showed 
that uniform background is not necessary. Our expla- 
nation needs only one assumption: f(£) vanishes at 
a certain finite energy. This assumption is quite nat- 
ural for the remnant of the merging of galaxies with 
MBHs, since MBHs would heat up the stars through 
the back reaction of the dynamical friction from the 
stars to the falling MBHs. 




Fig. 2. — The profiles of tlie binding energy dis- 
tribution N{£) of the galaxy after the BH settled 
in the center, obtained by our A^-body calculation 
( Nakano fc Makino 1999| ), where £ is the binding en- 
ergy, £ = ^ — ti^/2, and '^ is the potential energy 
defined to be minus the conventional gravitational 
potential. For all runs, the initial galaxy model is 
an isotropic King model with non-dimensional cen- 
tral potential Wc = 9. The mass of the BH, Mbh, is 
1/24 of the mass of the galaxy. 



3. Self-consistent Model of Weak Cusp 

Figure |k shows the self-consistent solutions of the 
density profiles of the galaxy with the central BH, 
for three different forms oi N{£); a King model with 
Wc = 9 (solid), a constant N{£) (dot-dashed) and 
exponential with cutoff N{£) ex (e^"^^ — 1) (dot- 
ted), where £o is cutoff energy. The mass of BH is 
l/24Mgai. We used the iterative method introduced 
by Binney (1982) to obtain these self-consistent solu- 
tions. Figure yb is the same result as Fig.ya but using 
the Hernquist model as initial guess of iteration. In 
all cases, the steep outer slope switches to the weak 
cusp with p ex r~^" at around r — 0.1. We can 
clearly see that the difference in the functional form 
of N{£) does not affect the slope of the central cusp. 
This result is in good agreement with our analytical 
explanation described above. 

Figure |j shows the solutions of the density pro- 
files with different BH masses. The size of the weak 
cusp region (or the so-called "break radius" ) is larger 
for larger BH mass, but the slope of the weak cusp 
remains unchanged. 

The relation between BH mass and size of weak 
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Fig. 3. — (a) Self-consistent solutions of the density 
profile fo r various N{ £), obtained by the iterative 
method ( Binney 1982 ). (b) Same results as (a) but 
using the Hernquist model as initial galaxy model. 



cusp region (hereafter simply term it "core") can be 
understood as follows. Suppose that the inner region 
of initial power-law density profile p = Kir^ is trans- 
formed by sinking BH to a shallow cusp p — K2r'' 
with radius r^ where —3 < p < q < 0. Strictly 
speaking, a King model has a small flat core, but we 
can neglect its contribution in the following discus- 
sion. Roughly speaking, the total mass of the region 
affected by the BH must be about the same as that 
of BH. Thus we have 



Mbh 



K2r'^ ■ 47rr dr : 



KirP ■A-rrr'^dr. (4) 



Using K2r^ = Kir^, we obtain 



(5) 



If the initial profile is isothermal {p = —2), the core 
radius would be proportional to BH mass. 

Such correlation between BH mass and core size is 
consistent with the observations of ellipticals, which 
indicate that there are linear correlations between the 



core radius and the effective radius ( Faber et al. 1997 ) 
and between the central BH mass and total mass 



dMagorrian et al. 19981) 
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Fig. 4. — Self-consistent density profiles with different 
BH masses. The initial model is the King model with 



4. Summary 



In this paper, we showed that the density cusp 
around MBH has the slope of -1/2 if the MBH has 
"fallen" to the center in dynamical timescale. Thus, 
we now understood why such weak cusps were formed 
i n numerical simulations of merging with central BHs 
( Makino fc Ebisuzaki 1996) o r simulation of sinking 
BHs ( iNakano fc Makino 19991 ). 

This is the only known process to form a weak 
cusp in the center of galaxies. Therefore, our result 
very strongly suggests that luminous elliptical galax- 
ies with weak cusps have experienced such "fall-in" of 
MBHs. As suggested by ME96, mergings of galaxies 
with MBHs already in the center is the most natural 
scenario for such an event. In other words, luminous 
ellipticals are most likely merger remnants. Our result 
also gives us an important suggestion for the origin 
of central MBHs, which are believed to exist in many 



galaxies (Richstone et al. 1998). The central MBHs in 
luminous elliptical galaxies with weak cusps were not 
formed there by some process such as gas accretion 
but were imported dynamically from the progenitor 
galaxies through the merging of them. If the cen- 
tral BH was formed in the timescale longer than the 
dynamical timescale, the central density cusp would 
have the slope of —3/2 or steeper (Bahcall & WoU 
1976; Young 198C|). Therefore, the observed shallow 



cusp and the existence of central MBHs are consis- 
tent only if ellipticals with shallow cusps are merger 
remnants. The clear dichotomy of weak cusps and 
steep cusps, and the correlation between the slope 
of the cusp and the normalized rotation velocity v/a 
(Figure ||) suggests that not only the merger of gas- 
poor ellipticals but also the "major merger" ([Barnes 



fe Hernquist 1992) of spirals (Barnes 1997) resulted 
mostly in weak cusps. Note that this clear separation 
is also visible in the sample of CaroUo et al (1997, 
Figure 8). This connection between central slope and 
kinematics implies that in the merging of two spirals, 
BHs would have become massive before two galax- 
ies finally merge, through gas-fueling induced by the 
tidal torque ( Selwood fc Moore 1999| ), and moreover, 
the rest of the gas would be ejected or formed into 
stars more rapidly than the BHs sink to the center of 
merger remnant. 

We thank Yoko Funato for useful advice on the nu- 
merical method and for valuable discussions on our 
results. We are grateful to Toshiyuki Fukushige for 
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Fig. 5. — Normalized rotation velocity of galaxies v/a 
versus cusp slope of luminosity densi ty profile at r = 
0.1" o btained by HST ob servations ( Gebhardt et al. 
19c(6|; iFaber et al. 1997| ). Filled circles denote the 



galaxies brighter than My — —21 and open circles 
denote those fainter than My = — 21. 
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